Non-Markovian stochastic Schrodinger equations: Generalization to real-valued noise 

using quantum measurement theory 



Jay Gambetta and H. M. Wiseman* 
Centre for Quantum Dynamics, School of Science, Griffith University, Brisbane Alll, Australia 

(Dated: February 1, 2008) 

Do stochastic Schrodinger equations, also known as unravelings, have a physical interpretation? 
In the Markovian limit, where the system on average obeys a master equation, the answer is yes. 
Markovian stochastic Schrodinger equations generate quantum trajectories for the system state con- 
ditioned on continuously monitoring the bath. For a given master equation, there are many different 
unravelings, corresponding to different sorts of measurement on the bath. In this paper we address 
the non-Markovian case, and in particular the sort of stochastic Schrodinger equation introduced by 
Strunz, Diosi, and Gisin [Phys. Rev. Lett. 82, 1801 (1999)]. Using a quantum measurement theory 
approach, we rederive their unraveling which involves complex-valued Gaussian noise. We also de- 
rive an unraveling involving real- valued Gaussian noise. We show that in the Markovian limit, these 
two unravelings correspond to heterodyne and homodyne detection respectively. Although we use 
quantum measurement theory to define these unravelings, we conclude that the stochastic evolution 
of the system state is not a true quantum trajectory, as the identity of the state through time is a 
fiction. 

PACS numbers: 03.65.Yz, 42.50.Lc, 03.65.Ta 



I. INTRODUCTION 

In nature, a quantum system is most likely found in 
an entanglement with at least one other quantum sys- 
tem. An example of this is a two level atom (TLA) im- 
mersed in an environment of harmonic oscillators (the 
electromagnetic field). This type of quantum system, a 
small system interacting with a larger system (the bath) 
is called an open quantum system [1]. The system-bath 
interaction causes the two systems to entangle, resulting 
in a combined state whose evolution can be theo- 

retically determined by the Schrodinger equation. How- 
ever, due to the many degrees of freedom of the bath, this 
is generally impractical and it is best to describe the sys- 
tem (TLA) by the reduced state p Ycd (t). The evolution 
of Pred(i) is found by averaging the outer product of the 
Schrodinger equation over all the possible bath states, 

ft«i(t) = TWI*(t)><*(t)|]. (i-i) 

Under the Born-Markov approximations [2] it is pos- 
sible to obtain a closed equation for p Iod (t). For mathe- 
matical consistency, this should be of the Lindblad form 
[3] . If there is a single Lindblad operator L (such as the 
lowering operator for the system) then this is an equation 
of the form 

A-cdW = -i[H,p Ied (t)]+jV[L]p Ted {t), (1.2) 
where H is the Hamiltonian and 

T>[L]p Icd = Lp rcd D - \DLp rcd - \p rcd DL. (1.3) 
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However this is only an approximation, in the non- 
Markovian situation in general one can not solve p red (t) 
or |\l/(t)) analytically, so p rcd (t) is difficult to determine. 

A breakthrough in solving this problem was achieved 
with the development of non-Markovian stochastic 
Schrodinger equations (SSEs). These stochastic dif- 
ferential equations for a state vector were first intro- 
duced for Markovian open quantum systems in mathe- 
matical physics [4, 5, 6, 7, 8, 9, 10, 11, 12] and then 
independently in quantum optics [1, 13, 14]. This ap- 
proach has subsequently been generalized to deal with 
non-Markovian systems [15, 16, 17, 18, 19]. In this pa- 
per we will follow the approach of Diosi, Strunz and Gisin 
(DSG) [19, 20, 21, 22]. In their approach, the system 
state vector \vj z (t)) [44] depends upon some (not neces- 
sarily white) noise z(t), which is drawn from some proba- 
bility distribution. The SSE has the property that when 
the outer product of \tp z (t)} is averaged over all the pos- 
sible z(t) one obtains p led (t). That is, 

p red (i)=E[|^(i))(^W|], (1.4) 
where E[...] denotes an ensemble average over all possible 

z(t)'B. 

In cases where an exact non-Markovian SSE can be 
derived, it is also possible to find an exact solution 
for p Icd (t) by other means. A key advantage of non- 
Markovian SSEs lies in the cases where no exact solution 
is possible. In this case approximations must be made 
in either approach. The advantage of the SSE approach 
is that the ensemble average Pi ed (t) is, by construction, 
guaranteed to be a positive operator. This fundamental 
property of a state matrix is not guaranteed by other ap- 
proximate equations for p Ied (t). This is true even in in 
the Markov limit; quantum Brownian motion is a case in 
point [23] . The other advantage of the SSE approach in 
general is that it allows the evolution of large systems to 
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be simulated numerically. This was the original motiva- 
tion for their introduction in quantum optics [13, 14]. 

Leaving aside the potential usefulness of SSEs, one may 
ask the question: is there a physical interpretation for the 
solution of an SSE, or is it simply a numerical tool for 
finding p TC d(t)? In the Markovian limit, that is when the 
master equation has the form Eq. (1.2), the answer is yes. 
The solution to the SSE, termed by Carmichacl is a quan- 
tum trajectory [1] , it can be interpreted as the state of the 
system conditioned on the measurement results obtained 
by continuously monitoring the bath [24] . For the Marko- 
vian case, different sorts of SSEs exist. They may involve 
jumps or diffusion, and are termed different unravelings 
of the master equation [1]. These different unravelings 
correspond to different detection schemes, such as pho- 
ton counting [1, 13, 14], homodyne [1, 24, 25], and hetero- 
dyne [25] detection. Other generalizations [26, 27, 28, 29] 
have also been investigated. 

In this paper we will investigate the question of phys- 
ical interpretation of non-Markovian diffusive SSEs of 
Diosi, Strunz and Gisin (DSG) [16, 19, 20, 21, 22]. 
We will show that quantum measurement theory (QMT) 
does give meaning to the | - z (t)) at any particular time, t. 
However, the linking of the state [^(i)) at different times 
to make a trajectory appears to be a convenient fiction. 
We also show that the theory of DSG can be general- 
ized by considering different sorts of measurements (un- 
ravelings) on the bath. We use our approach to define 
two different unravelings. The first results in DSG SSEs, 
with complex- valued noise z(t). In the Markovian limit 
this unraveling corresponds to heterodyne detection. The 
second, which can only be defined for some system-bath 
couplings, has real- valued noise and has homodyne de- 
tection as its Markovian limit. 



II. SYSTEM DYNAMICS AND QUANTUM 
MEASUREMENT THEORY. 

A. Schrodinger Equation for the combined system 

With h = 1, a system interacting with a reservoir of 
harmonic oscillators has the total Hamiltonian 

H tot = H Q + H + ff bath + V. (2.1) 

Here the system Hamiltonian has been split into Hq (the 
action of which is described later) and H (the remainder) . 
The Hamiltonian for the bath is 

-ffbath = y^fcajflfc, (2.2) 

fe 

where k labels the modes of the bath, dk and tOk are the 
lowering operator and angular frequency of the k th mode 
respectively. We assume the interaction Hamiltonian to 
have the form 

V = i(W - btt), (2.3) 



where L is a system operator and where we have define 
the bath lowering operators b as b = J2k9k a k- That is, 
the coupling amplitude of the k th mode to the system is 
9k- 

The Schrodinger equation for the combined state is 

dt\V(t))=-iH tot \y(t)), (2.4) 

which can equivalcntly be written as 

|*(t))=t/(i,0)|*(0)>, (2.5) 

where U(t, 0) is called the unitary evolution operator. 
Defining a unitary evolution operator for the 'free' system 
and bath as 

U (t, 0) = e - l ^ 0+i ^ h )(*-°>. (2.6) 

We can write U(t, 0) as U(t, 0) = U (t, 0)U int (t, 0), where 
Uint{t, 0) is the unitary evolution operator that describes 
the total evolution with the free dynamics removed. 
We can then define an interaction picture state as 

|*int(*)> =EWt,0)|*(0)>, (2.7) 

which obeys 

dt|*«it(*)) = -i(H iat (t) + V^ntW) l*intW)- (2-8) 
The Hamiltonians in the interaction picture are 

Hint(*) = U^t,0)HU Q (t,0), (2.9) 

and 

v int (t) = i(SL(*)£tat(*)-Sint(*)iL(t)). (2-io) 

where 

km(t) = Y,9ka k e- iUkt , (2.11) 

k 

L int (t) = Le~ 1 ^ 1 . (2.12) 

Here we have finally restricted Hq to be such that L in the 
interaction picture simply rotates in the complex plane 
as indicated in Eq. (2.12). The interaction picture can be 
viewed as moving the time dependencies due to the free 
bath and system dynamics from the state to the opera- 
tors. Unless otherwise stated the rest of this paper will 
be in the interaction picture and thus we will drop the 
subscripts 'int'. 

B. QMT and Conditional System States 

In open quantum systems a measurement is always per- 
form on the bath. Due to the entanglement between the 
bath and the system the measurement on the bath re- 
sults in an indirect measurement of the system [30] . The 
state of the system after the measurement is dependent 
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on the results of the measurement, so we call this a con- 
ditional system state. To mathematically describe this 
(for a more detailed description see [26, 30, 31]) we de- 
fine |{(?fe}) as the arbitrary basis the measurement is per- 
formed in. Note that |{<?fc}) does not necessarily have to 
be normalized. For our purposes we restrict |{(7fc}) to be a 
state in the interaction picture with no time dependence 
(it will be Uq (t, 0)\{qk}) hi the Schrodinger picture). A 
typically example of this is a coherent bath state. This is 
the state (in the interaction picture) the bath (harmonic 
oscillators) has when driven by a classical current [32]. 

In the basis we can define a probability- 

operator-measure (POM) element, or effect, as 



F{ Qk} = \Uk})({q k }\- 



(2.13) 



Here the subscript {qk} is the result of the measurement. 
The effect is important as it allows one to calculate the 
probability density of results {qk} 



P({q k },t) = (*(t)\F {qk} \*(t)}. 



(2.14) 



If one were only interested in obtaining probabilities the 
effect would be all one would need. However, since we are 
interested in the state of the system after the measure- 
ment, we need to define a set of measurement operators. 
The constraint the measurement operators must obey is 
rt 



F{ qk y = M^ qk yM{ qk }. For example, we can decompose 
the measurement operators as 



M. 



{q k } = \i n k}) ({qk}\, 



(2.15) 



where {n k } is arbitrary, and is the state the bath is left 
in after the measurement. Since in most detection situa- 
tions a measurement results in annihilating the detected 
field the most natural choice for {rife} is the vacuum state 

{04- 

In QMT the combined state after a measurement at 
time t, which yielded results {qk} is [30, 31] 



!*{**>(*)) 



M {gk} Mt)) 



(2.16) 



Using equation (2.15), with n k = for all k, the 
combined state after the measurement is, \^[ qk y(t)} = 
|{0fc})|^ {9fc }W), where 



(2.17) 



Equation (2.17) is the conditional system state and we see 
here directly how the entanglement between the bath and 
the system results in the system state collapsing upon 
measurement of the bath. One of the properties of this 
conditional system state is that p le d(t) (equation (1.1)) 
can be written as 

AedW = [({q k }Mt))(*(t)\{qk})d{q k } 



P({q k },tM {qk} (t))(iP {qk} (t)\d{q k } 

m{ qk }(t))(*P {qk} (t)\i (2.i8) 



where E denote an average over the distribution 
P({q k },t). From equation (1.4) we see that the con- 
ditional state satisfies the same requirements as a so- 
lution of a SSE. This suggests that the time derivative 
of equation (2.17), if it could be written in terms of 
l^{<?fc}W)' cou ld be interpreted as a SSE. One problem 
in determining the time-derivative is that Eq. (2.17) in- 
volves the probability P({q k }, t), which requires knowing 
IvE'(t)), and, as mention earlier, this in general is indeter- 
minable. However, this problem may be overcome using 
linear quantum measurement theory (LQMT). 

LQMT uses the same principles as QMT except we use 
an ostensible distribution (A({qk})) in place of the actual 
probability [26, 33]. As its name suggests, the ostensi- 
ble probability distribution need bear no relation to the 
actual probability distribution. However, it must be a 
proper probability distribution (non-negative, and inte- 
grating to unity), and must be non-zero wherever the ac- 
tual distribution is non-zero. Using the ostensible prob- 
ability distribution, the conditioned system state is 

fc)W) = iiM . (2 , 9) 
v A (fel) 

We will call it the linear conditioned system state, be- 
cause it depends linearly on the pre-measurement state 
|*(i)), unlike Eq. (2.17). Since A({q k }) is not equal 
to the actual probability, \ip{ qk }(t)) will not be normal- 
ized and to signify this we use a tilde above the state. 
Note that this notation, following our earlier convention 
[26, 28], is the reverse of that used by DSG [21]. Because 
it is unnormalized, the linear conditioned system state 
does not have a clear physical interpretation. However, 
it still is useful as it is easier to calculate (involving only 
linear equations), and we can write 



Pred(t) 



({q k }mt))m)\Uk})d{q k } 

M{qk})W{ qk} (t))ti {qk} (t)\d{qk} 



E[fc } (t))«; { , fc }(*)i] 



(2.20) 



where E denote an average using the ostensible distribu- 
tion A({q k })- The condition for obtaining a linear SSE is 
we have to be able to write the time derivative of equa- 
tion (2.19) in terms of only \ip{ qk y(t)}. 

A linear SSE is only really useful if it can be trans- 
formed into a nonlinear SSE for the normalized state 
\4 } {q k }(t)) • To do this one requires that there exists a 
Girsanov transformation for the variables {qk} [34]. This 
is a transformation that takes into account the relation 
between the actual probability to ostensible probability, 

P({q k },t) = $ {qk} (t)\iP {qk} (t))A({q k }), (2.21) 

which follows from Eqs. (2.19) and (2.14). Specifically, 
the Girsanov transformation is a time-dependent trans- 
formation that changes the variables {qk} into the vari- 
ables {q£} such that 



A({<z£})dK A } = P(tik},t)d{q k }. 



(2.22) 
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We can see the usefulness of this transformation as 
follows. If we normalize the unnormalized states, but 
keep the same ostensible distribution, then the ensemble 
average will not reproduce p le d(t)- 

^^\^ {qk} (t)){i> {qk} (t)\d{q k } ± p Ied (t). (2.23) 

However if {q k } are chosen from the actual distribution 
then, of course it does: 



n l qk ^\ ^ {qk }(t)){^ qk }mM = PrcdW- (2.24) 



Equivalently, using the ostensible distribution for {q^}, 



A ({^}),7 



I^Wl 



l^{,»}(*)>^{«»}(*)l«} = Pred(i). (2.25) 



Note that both {q k } and {q£} appear here. This means 
that if we have a linear SSE, we can derive a nonlinear 
('actual') SSE by normalizing the state 



IVw(*)> 



1 



where 



(2.26) 



(2.27) 



but generating the SSE by drawing {q£} rather than {qk} 
from the ostensible distribution. 

Now that we know how to use Eq. (2.26), we can 
calculate the time derivative of \ipt qk }(t)) in terms of 

l^{ «*}(*))• This res ults in 

dtW {qk} (t)) = -J—d t \i> {qk} (t)) + \^ {qk} (t))d t -l- 

(2.28) 

where 

d t\^{ qk }(t)) = dt\i>{q h }(t)) +^2d t qkd qk \^ {qk} (t)). 

k 

(2.29) 

Here we have assumed that we can define dtqk so as to 
generate a qk(t) which ensures that Eq. (2.22) is always 
satisfied. From the above discussion, it is thus apparent 
that three conditions must be satisfied if Eq. (2.28) is to 
be a SSE for the system state \tl>{q k }(t)). These are: 

1. It is possible to obtain a linear SSE. that is 

W{»}(*)>- 

2. There is a Girsanov transformation {q£ } — > {qk(t)} 
such that an equation for dtqk for all k can be found 
explicitly 

3. Equation (2.28) can be written in terms of only 

IVw(*)>- 

If we can satisfy all these conditions then we have a 
SSE which generates a state with a definite physical in- 
terpretation. The SSE generates a state at time t which 



is of the form of Eq. (2.17). This is clearly the normalized 
state conditioned on a measurement being performed at 
time t on the entire bath, and yielding results {qk}- 

It is important to note, however, that the linking of the 
states at earlier times to form a trajectory (which is how 
the SSE generates the state at time t) appears to be a 
convenient fiction. A measurement on the whole bath at 
time t is clearly incompatible with a similar measurement 
at an earlier time. It is only in the Markovian limit that 
compatible bath measurements can be made, so that the 
quantum trajectory as a whole can be interpreted phys- 
ically. In other words the time evolution generated by 
the SSE simply links together hypothetical conditioned 
states at different times, with different measurement re- 
sults {qk(t)}. The relation between the results at dif- 
ferent times is purely mathematical, not physical. The 
mathematical relation comes from the time-dependent 
Girsanov transformation: the q£ corresponding to the 
qk(t) are the same at all times. 



III. COHERENT BATH UNRAVELING 
A. Coherent Noise Operator 

The first unraveling we consider is that associated with 
the bath being projected into a multitude mode coherent 
states, that is |{sfe}) = |{&/c}) where 




Note that these states are deliberately not normalized, 
so that the multi-mode integral of the effect -F!r 0fc } = 
K a fe})({ a fc}| is unity. We call the resultant unraveling 
the 'coherent state unraveling'. For this unraveling we 
define the noise operator 



E 



SfcOfce 



(3.2) 



where O/. = ujk— u>q. This noise operator has the property 
z(t)\{a k }) = z(t)\{a k }}, (3.3) 
where z(t) is the noise function, given by 

z(t) =J29ka k e- mkt . (3.4) 

k 



An important property of the bath is its correlation: 
how the noise operator (function) at time t is related to 
that at time s. This is determined by the commutator 
(operators) or correlation function (noise functions). For 
a non-Hcrmitian operator there are two important com- 
mutators, 



[*(*).*(*)] 
[z(t),z( S )t] 



0. 

a(t — s), 



(3.5) 
(3.6) 
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where, in the notation of DSG, 



a(t — s) 



E 

k 



9k\ 2 e 



-iQ.k{t—s) 



(3.7) 



which we call the memory function. 

The second form of correlation in defined in terms of 
the noise functions as E[z(t)z*(s)]. This depends on the 
probability for obtaining the results {a k } in the measure- 
ment at the two times. In linear QMT, these probabilities 
are given by the ostensible distribution A({a k }), which 
may be chosen to be time- independent. It is convenient 
to choose A({cifc}) to be equal to the actual probability 
that would arise when the bath is always in the vacuum 
state. That is, 

A({a fe }) = ({0 fe }|H}><H}|{0 fc }) =7r -« e -£>*l 2 , 

(3.8) 

where k — ^2 k - As will be seen later, this is appropriate 
if the bath is initially in this state. The correlation for 
the noise functions under this assumption is, 



E[z(t)z*(s)] 
E[z(t)z(s)} 



0. 



(3.9) 
(3.10) 



Note that we have used the notation discussed below 
Eq. (2.20). Thus for the special case where the ostensible 
probability is given by Eq. (3.8), the memory function is 
equal to the correlation of the noise functions. 



1. The Markov Limit 

Since one of our aims is to consider the Markovian 
limit of our non-Markovian SSEs (in which one obtains 
a genuine quantum trajectory), the Markov limit of all 
our main results will be presented. In the Markov limit 
the number of modes become continuous and the cou- 
pling constant \g k \ becomes fiat (|^| = g) and equal to 
^ j/2n . This allows us to write 

7 / ^-i( u - UQ )(t- s ) ^ 



a(t — s) 



2^ jo 
j_ 
2tt 



(3.11) 



and for optical situations (high uiq situations) with little 
error this can be written as 



a(t — s) = 



2tt 



-^-^dn = 7<y(t- s). (3.12) 



Therefore, 

E[z(t)z*(s)} = [z(t),z( s y] = jS(t-s), (3.13) 
E[z{t)z(s)] = im,z(s)} = 0. (3.14) 

This implies that ostensibly z(t) is a complex gaussian 
random variable (GRV) of mean and variance "f/dt. 
That is, z{t) = ^((t), where ((t) is the standard com- 
plex white noise function [35] . These are the correct cor- 
relation function for the heterodyne noise functions [26]. 



B. The Linear Stochastic Schrodinger Equations 
for the Coherent Unraveling 

In this section we will derive the linear non-Markovian 
SSEs for the ostensible probability introduced above, and 
show that in the Markov limits it gives the linear Hetero- 
dyne SSE. We use many of the same techniques as DSG. 
To calculate the linear SSE we write the Schrodinger 
equation in terms of the noise operator, z(t) 



|*(t)) = {-iH(t) + z\t)L - $(t)#}|tf(t)>. (3.15) 



Then by differentiating equation (2.19) with respect to 
time (with q k set to afc) we obtain 



9t|Vw(*)) = {-iH{t) + z*{t)L)$ {ak) {t)) 

VMM) ' 



(3.16) 



as H(t) is a system-only operator and ({afe}| is the left- 
eigenstate of z(t)^. To satisfy the condition for a linear 
SSE we must evaluate the last term in this equation in 
terms of \i/j{ ak y(t)). To do this we use [32], 

{{a k }\a k Mt)) = (^ + 5„;)({a Jfc }|*(t)) (3.17) 

and 

da'({a k }\^(t)) a k - 
V A (l a fel) z 
With these two expressions and the definition of z(t), 

v A li a *;) k 

This allows us to write equation (3.16) as 



dS{a k] (t)) = {-iH(t) + «*(*)L-Lt^ fffc e-* 



— iQut 



:0-;}l^w(*)>. 



(3.20) 



This is a linear equation in terms of {a k }- Note that 
it is not really a SSE, as the final term implies that the 
evolution of the state \i/j{a k }(t)) depends not only on it- 
self, but upon neighbouring states with different values 
of {a k }. That is, we cannot simply choose (stochasti- 
cally) a value for {a k } from the ostensible distribution 
and then propagate forward the system state using that 
value. However, we can make progress towards an equa- 
tion where we can do this by rewriting the partial deriva- 
tive in terms of a functional derivative. This is done by 
using the following relation (see for example [37]), 



S dz*(s) 



5z*(s) 



ds , 



(3.21) 
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where is the initial time. This gives 

d t \Mt)) = [-iH{t) + z*{t)L-D J a{t-s) 

x s^) ds }^* {t)} > (3 ' 22) 

where a(t — s) is defined in equation (3.7). By replac- 
ing the partial derivatives by the functional derivative we 
have enforced the initial condition |W(0)) = |{0fc})|?/ , (0)), 
This is seen as follows. At t = the functional deriva- 
tive term in the above equation will have zero contri- 
bution, from the definition (3.2f). By comparison with 
the corresponding term in Eq. (3.20), it follows that 
9al\"<P{ ak }(t)}\t=0 = for all fc. From Eq. (2.19) this 
is only possible if the system and bath states initially 
(at time 0) factorize, and if A({a fc }) = | ({a fc }|V>bath) | 2 - 
From our choice (3.8) of ostensible probability, this en- 
forces IV'bath) = |{0fe})- This is physically acceptable as 
we may assume that at time the system and bath are 
uncoupled, and the bath is in the vacuum state. 

Like Eq. (3.20), Eq. (3.22) is not really a SSE be- 
cause the functional derivative means that it depends 
not upon a state \ip z (t)} at all times for a single value of 
the function z(t), but rather also upon states for other 
values of that function. That is, we cannot stochastically 
choose z{t) in order to generate a trajectory independent 
of other trajectories. Instead, all possible trajectories 
would have to be calculated in parallel. This means that 
the amount of calculation involved in solving Eq. (3.22) 
would be comparable to that required for directly solving 
the Schrodingcr equation (2.4). However in some circum- 
stance we can make the following Ansatz [21], 

^_|^))=4(f, s )|^(t)}, (3.23) 

where O z (t,s) is some system operator which is a func- 
tion of t, and s, and a functional of z. With this Ansatz 
the linear SSE becomes 

fit $*(*)) = {-iH{t) + z*{t)L-D j a(t-s) 

xd z {t,s)ds}\Mt)). (3.24) 

This is now a true SSE, where each trajectory can be 
evolved independently. It is the same as the linear SSE 
that DSG presented in Ref. [20, 21]. Note that it is 
non-Markovian because the noise z*(t) is non- white, be- 
cause of the finite lower limit of the integral, and because 
O z {t, s) may depend upon z. 



1. The Markov Limit 

The next question is what is the Markov limit of this 
equation? To find this we use the results of section (III A) 



and the fact that O z (t,t) = L [21]. Applying them to 
equation (3.24) results in 

dt\Mt)) = {-iH(t) + z*(t)L - T £tLj|^( t )) ) (3.25) 

where z(t) = y/j((t). By its method of derivation, this 
equation is in Stratonovich form [35]. To compare with 
the standard Markov equations we should convert it to 
an Ito SSE. This can be derived by using an arbitrary 
basis and defining i/jj = (j\ip) and Lj t k = (j|L|fc). Then 
if the Stratonovich form is 

d t ipj = cij +b j (*(t), (3.26) 

the Ito form (which we indicate by use of the infinitesi- 
mals rather than the derivatives) is 

dt & 

= a dt + b,dC(t)dt + ^J2 h *^ b r ( 3 ' 2? ) 

The final term here is the Ito correction term. Looking 
at equation (3.25) we see that, bj = Ylk Lj^k, an d 
since dipk/dij)* is zero for all fc, the correction term for 
this equation is 0. Thus the Ito SSE is, 

d\v} z (t)) = dt(-iH(t)+Lz*(t)-^pL)$ z (t)), (3.28) 

which is the standard linear heterodyne SSE presented 
in Ref. [33] as z(t) = ^ Q{t) = V'y (&(*) + *'&(*)), where 
£fc(i) are the standard real- valued white noise terms [35]. 



C. The Actual Stochastic Schrodinger Equations 
for the Coherent Unraveling 

In this section we will derive the non-Markovian SSEs 
for the actual probability distribution and show that in 
the Markov limits it gives the the usual heterodyne SSE. 
Again, we use many of the same techniques as DSG. 

As discussed in section II B, to find an actual (i.e. non- 
linear) SSE for the normalized state we need to satisfy 3 
conditions. The first was to derive a linear SSE, which we 
did in the preceding section (by making use of an anstaz). 
The second condition is to find random variables with 
the actual probabilities of measurement results. To work 
out these random variables, {a^} we use the Girsanov 
transform (2.21) to find a first-order partial differential 
equation (PDE) for the probability, from which the char- 
acteristic equation generates the transformed variables. 

To obtain the PDE we differentiate Eq. (2.21), giving 

d t P({a k },t) = (Tjj {ak} (t)\d t \yl {ak} (t))A({a k })+c.c. 

(3.29) 

By equation (3.20) the above becomes 

d t p{{a k },t) = {(v; W} (t)|L|v; W} (t))^4^ e ^ i 

fc 

+c.c.}A({a fc }), (3.30) 
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Using the fact that \fp{a k }(t)) is analytical in a* k (so that 

d ak IV>{a fe }(*)) = 0) [20], and the product rule for differ- 
entiation, we can simplify the above to 

d t P{{a k },t) = -J29ke- mkt d al {$ {ak} (t)\V 

k 

x|$ {0fc} (t))A({o fc })} + c.c. (3.31) 

Defining 



(^)t = (<P {ah} (t)\L^ {ak} (t)) = 



^ {ak} {t)\L^ {ak] (t)) 

<^W}(*)IVW(*)) 

(3.32) 

allows us to write 
d t P({a k },t) = -^ 5fe e-^^ a .{(it) t P({ afe}) t)} 

(3.33) 



k 

+c.c. 



This is the PDE for the probability distribution. 
At t = 0, we have from Eq. (2.21) that 

P({a k }, 0) = $ {ak} (0)\i> {ak} (0))A({a k }). (3.34) 

As noted above, to obtain Eq. (3.22) we had to assume 
that the bath was initially in the vacuum state, uncor- 
rected with the system. This enforces the equation of 
the initial probability distribution to be the ostensible 
distribution 

P({a k },0) = A({a k }) =^- K e -E,.KI 2 . (3.35) 
From this PDE we can find the characteristic equations 
d t al=g k e- in * t (ti) t , (3.36) 



which integrates to give 



4(*) = 4(o) 



(3.37) 



The random variable a£(0) is one with probability distri- 
bution (3.35). With equation (3.37) and our noise func- 
tion definition, Eq. (3.4), we can write z(t) as 

z*{t) = 4(0)^e^* + f a*(t - s)(fr} s ds. (3.38) 



The term a£(0)^e lOfc * is the noise function one would ob- 
tain if the bath were assumed to be in the vacuum state. 
This is our assumption for the ostensible distribution so 
we will label this term z*^(t). This allows us to write 

z*(t)=zl(t)+ [ a*(t-s)(D) s ds, (3.39) 
Jo 

where z^(t) obeys the correlations expressed in equations 
(3.9) and (3.10). 

The third condition was to show that we can write 
equation (2.28) in terms of only \ip z (t)). To do this we 
start by calculating dt\ip z (t)). Using equations (2.29), 
(3.22) and (3.21) we get, 



dSz(t)) = {-iH(t)+z*(t)L-{P-(fr) t ) 

-ds}\i> z (t)}. (3.40) 



; / a(t — s 
lo 



S 



Looking at Eq. (2.28) we see that to obtain the actual 

SSE we need to calculate \i>z(t))dt\'ip{a k }\ ■ Using the 
above, 



\Mt))dt 



\^{a k }\ 



\Mt)) 

\^{a k }\ 



«Vz(i)K|^(i)> 



{z*W(L) t -(^(t)|(it_(Lt )t) _J_ 



+C.C 



Therefore Eq. (2.28) becomes 
dt 



\Mt)) = {-iH(t) + z*(t)L}\Mt))-(V-(V) t )-J— f t a(t- S ) T 4 r ^ 

1 J W{a k }\J0 



8z*{sy 



\Mt)){**V)(L}t - ®,®\<P ^yi~ 

}/ 



-c.c /2 



(3.41) 



(3.42) 



This can be simplified by using the fact that if our SSE has the form dt\ip) — (A + B/2 + B*/2)\ip) then we can 
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define a state \(j>) = exp(/ (B — B*)dt/2)\ip) (which is the 
same state as that gives a equivalent SSE, of form 



J 



(A + B)\<fi). Applying this to the above gives 



t \Mt)) = {-*H(t) + z^t)(L-(L) t )\\Mt))-(V-(V) t )-J— fa{t-s) 

1 J \4>{a k }\ JO 



Sz*{s) 



\^{a k }\ Jo 



t g 



(3.43) 



r 



This is not yet a SSE as it still contains \ip z (t)) 
terms, however if we can make the Ansatz described by 
Eq. (3.23) we can write this as 



dt\Mt)) 



iH(t)+z*(t)(L-{L) t ) 

a(t -*){(# -(fr) t )6 z (t,s) 
(V-{tt) t )6 z (t,8))}d8\\4>,(t))> 



(3.44) 

which is a genuine SSE. This means that an actual SSE 
(generating normalized states with their actual proba- 
bilities) can only be found if we can make the Ansatz 
describe in Eq. (3.23). 

This SSE is the same as that presented in Refs. [21, 22]. 
As shown here, it gives us the state the system would be 
in if at time t we performed a measurement in the coher- 
ent basis, and the result was z(t) as defined in Eq. (3.39). 
Note that this means that the result z(t) depends upon 
the system state at earlier times in the trajectory gener- 
ated by the above SSE. We have argued above that this 
linking of states at different times is a convenient fiction, 
but we see here that it is mathematically necessary in 
order to generate measurement results for a particular 
time with the actual probability. 



1. The Markov Limit 

Finally, we are again interested in the Markov limit of 
this SSE. Taking the Markov limit of the noise function, 
one obtains 



z*(t) = z* A (t) + l(L% 



(3.45) 



where = V7C*(*)- 

To apply the Markov limit to Eq. (3.44) we use 
a(t — s) — > 7<5(< — s) and O z (t, t) = L, resulting in 



dt\Mt)) = { 



J^ + (L-(L) t )(z*(t) + l(V) t ) 
<LtL) t )}|^(f)>, (3.46) 



which is in Stratonovich form. To convert this to an 
Ito SSE we have to calculate the Ito correction term in 
Eq. (3.27). For this equation, the correction term is 

which with Eq. (3.46) results in, 



{L^L) t + (D)t(L)t) \ip z {t)}, 



(3.47) 



d t \Mt)) = {-^ + (L-(L) t )z*(t) 



7 



{DL-L{D) t ))\4, z (t)). (3.48) 



This is the Ito SSE for the actual measurement prob- 
abilities. When we substitute in z*(t) from Eq. (3.45) 
we get the same heterodyne SSE as that presented in 
Ref. [12, 25]. 

Readers familiar with quantum trajectory theory for 
heterodyne detection may be puzzled by the factor of 
1/2 multiplying the deterministic contribution to z(t). 
This function is, according to the above theory, the result 
of measuring the bath at time t in the coherent state 
basis. But in the usual quantum trajectory theory [25] 
the measured (complex) heterodyne current at time t is 

I(t) = ^fC(t)+l(L)t. (3.49) 

which lacks the 1/2. Where does this discrepancy come 
from? To answer this question we have to consider the 
definition of a measurement, and in particular the time of 
the measurement. In quantum trajectory theory we must 
consider the measurement which conditions the state at 
time t as actually occurring at a time t + dt [26]. That is, 
the ^-correlated bath must be given a chance to interact 
with the system before the measurement is made. By 
contrast, in the above theory the measurement occurs 
exactly at time t. For a non-Markovian bath (with a 
finite correlation time) the difference between t and t + 
dt is infinitesimal. However in the Markov limit, this 
infinitesimal difference in measurement time causes the 
finite difference between z{t) and I(t). 

It is easiest to see this using the Heisenberg picture. 
From the above theory, 

E[*(*)] - (*(t)|i(t)|*(f)) 

= (m\({o k }\uUmt)Uint(t)\{Ok})\m) 
= (m\({0k}\z H (tw k })\m), (3.50) 
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where zji{t) is the Heisenberg noise operator. In quan- 
tum trajectory theory the measurement is defined to take 
place after the system and bath have interacted for a time 
dt, so that 

E[J(t)] = (y(t + dt)\z(t)\V(t + dt)} 

= (V>(0)|({0 fe }|c4(* + dt)z{t)Uint(t + dt) 

x\{o k })\m) 

= (^0)|<{0 fe }|/(i)|{0fc}>|V<(0)>. (3.51) 

Therefore, 

I(t) = U? at (t + dt,t)z H (t)U iat (t + dt,t). (3.52) 

By using standard Heisenberg equations it can be 
shown that 

/t+dt 
alt - s)Ul t {s)LU- mt {s)ds, (3.53) 

which has a Markov limit of the form 

7(f) = z H (t) + ^uUt)LU int (t). (3.54) 

This is the operator form of the Heterodyne current, and 
shows the extra contribution discussed above. It is simi- 
larly easy to show that the Markov form of zn{t) is 

*A(t) + \ul t {t)LU int {t), (3.55) 

where z^(t) = J^k 9k&k{0)e~ inkt . These relations are 
analogous to the Markovian input-output theory of Gar- 
diner and Collett [38] . The correspondences are as follows 

z A {t) ~ b in (t), (3.56) 
z H (t) <-> b(t), (3.57) 
/(*) b out (t). (3.58) 

IV. QUADRATURE BATH UNRAVELING 

In this section we will present a second unraveling 
which is conditioned on real noise and has homodync 
detection as its Markov limit. 



A. Quadrature Noise Operator 

To obtain a SSE with real noise, it is natural to con- 
sider a quadrature noise operator, 

z(t) = b{t)e lUot e- l,t > + &t(t) e -*"°* c * (4.1) 

where b(t) is defined in equation (2.11) and </> is some 
arbitrary phase. The noise operator has a two-time com- 
mutator 

[z(t),z(s)} = a(t - s) - a*(t - s), (4.2) 



independent of 4>. The phase <fi defines the measured 
quadrature: an x-quadrature measurement occurs when 
<f> is set to zero, and the conjugate measurement of the 
y-quadrature occurs when <f> = tt/2. Unless otherwise 
stated we will set <f) to zero. 

The basis for the bath measurement is |{gfc}) and must 
satisfy 

z(t)\{q k })=z(t)\{q k }). (4.3) 

The problem with this noise function is that it is hard 
(maybe impossible) to work out a time-independent 
eigenstate | {<?&}) in the interaction picture. However, 
we can find the eigenstate if we make the assumptions 
that for every mode k there exists another mode, which 
we can label —k, such that £l- k = — Cl k and g_fc = g k . 
These assumptions simply mean that the modes coupled 
to the system come in symmetric pairs about the system 
frequency loq. Without loss of generality we can take the 
<7fc's to be real, absorbing any phases in the definitions 
of the bath operators. With all of these assumptions we 
can rewrite equation (4.1) as 

z(t) = J2 2 9k (x+ cos(n k t) + Y~ sin(fM)) . (4.4) 

Here we have introduced the two-mode quadrature oper- 
ators 

X± = (x k ±x_ k )/V2, (4.5) 

= (y k ±y-k)/V2, (4.6) 
where x k and y k are the quadratures of d k : 

a k = (x k +iy k )/V2. (4.7) 
These operators have the commutators 

[X k ,Y k -]=i, [X k ,Y+} = 0, (4.8) 
[Xk,Yk}=0, [X+,Y+]=i. (4.9) 

Since {X^} and {Y k } form two mutually commuting 
sets of commuting operators, and thus have a common set 
of eigenstates. Since z(t) is a linear combination of these 
operators, the eigenstates of {X k } and {Y k ~} are the 
|{<?fc}) we seek. Therefore we can write the two eigenvalue 
equations, 

Y k \{q k }) = Y k -\{q k }), (4.10) 
X+\{q k }) = X+\{q k }). (4.11) 

This suggest that we should write | {<?&}) as Y^}), 
but for brevity we will continue to write it as |{5fc}). The 
form of the state that satisfies these equations, in the y k - 
basis, for a particular k is 

/ -^\{y' Y k -)/V2)- k \(y' + Y k ) / >y/2) k e~ iX i V, 
J v Zn 

(4.12) 
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while in the x^-basis it is 
dx 



/2tt 



|(X+ - x')/V2)^ k \(X+ + x')/V2) k e^ 



(4.13) 

Under these assumption we can show that the memory 
function a(t — s) in Eq. (3.7) becomes equal to the real 
function (3(t — s) given by 



p(t - s) = 2J2 \9k\ 2 cos(O fe (i- a)). 



(4.14) 



k>0 



Thus the commutator expressed in Eq. (4.2) becomes, 



[z(t),z(s)]=p(t-s)-P(t-s) 
Moreover, the noise function is 



0. 



(4.15) 



(t) = J2 2 9k (X + cos{£l k t) + Y k sm(Q k t)) . (4.16) 



fe>0 

Since X k and Y k are real, z(t) is also. 

We can define the correlation function for the noise 
functions as E[z(i)z(s)], and again this depends on the 
probability distribution for the variables X k and Y k . It 
is again convenient to choose the ostensible distribution 
to be that corresponding to the bath being in the vacuum 
state. Explicitly we then have 



A({X k ,Y k }) =tt 



= 7r -/2 P -E k>0 ( x t 2 +Y k - 2 ) 



(4.17) 



With the usual ostensible distribution the correlation 
function is 



E[z(t)z(s)\ = 2^2 \9k\ 2 cos(n k (t - a)) = /3(t - a) 



k>0 



while E[z(i)] = as before. 



(4.18) 



1. The Markov Limit 

The symmetry assumptions we have made in order to 
obtain this z(t) are compatible with the Markov limit 



in which the modes become continuous and the coupling 
constant becomes flat in fc-space (which of cause is sym- 
metric around ojo). As in the coherent case, the memory 
function (3(t — s) in the Markov limit equals "/5(t — s). 
Therefore in this limit the noise function is ostensibly 
given by z(t) = ^/y£,(t) where £(i) is a real- valued Gaus- 
sian white noise term [35]. 



B. The Linear Stochastic Schrodinger Equation for 
the Quadrature Unraveling 

To find the linear non-Markovian SSE we start by ap- 
plying our assumptions to the Schrodinger equation for 
the combined state 



dt|*(f)> = {-iH(t) + Y,g k [L(&l 



a_,,e 



') 



k>0 



-Zt^e-'^ 4 +a_ fe e^ t )]} | *(*)). (4.19) 
Now by Eq. (4.4) we can write this as 

dt|*(t)> = {-iH(t) + Lz-J29kL x {& k e- inkt 

k>0 

+e m > t a- k )}\*(t)}, (4.20) 

where L x = (L + ]J). Using definitions (4.5), (4.6) and 
(4.7) we rewrite the above equation as 

dt\V{t)) = {-iH(t) + zL-J29kL x (x+cos(n k t) 



fe>0 



+iY k cos(fifci) — iX k sin(Slfcf) 



-Y k sin 



in(n fc t))}|*(t)). 



(4.21) 



As in the coherent case to find a linear SSE we dif- 
ferentiate Eq. (2.19) with respect to time, except that 
this time |{(?fe}) is given by Eq. (4.13) and the ostensible 
probability is given by Eq. (4.17). Using Eq. (4.21) we 
obtain 



= (-ig(f) + *(*)£) fe>(*)) -^>£,{cos(^t)(z {{qkW ^ {t)) +X+\i> {qk} (t))) 

k>o ^JA({X+ 7 Y k }) 

+ ^4tfc } W)-»H)}- (4.22) 

I 



The inner products in the above equation can be simpli- fied to 



({ qk }\X k \^(t)) = iJLdq^t)), (4.23) 
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(tik}\Y k + \y(t)) = _i_({ 9fc }|*(f)). (4.24) 



as X k and Y^ have the commutators listed in equations 



(4.8) and (4.9). 
It can also be shown that 



d 



dYZ 



fc}(*)> 



d 



lV> te }(*)> = 



A({X+,Y k ~}) dYr 

+n"i^»}(*)>, 

1 d 



■<{?*}!*(*)> 



(4.25) 



A({x+,y,-}) ™* 

+*jffe}(*)>, 



+ <{«fc}|*(*)> 



(4.26) 



and using equations (4.23) and (4.24) with the above two 
equations we can write the inner products in terms of 
their conjugate variables. This allows us to write the 
linear equation as 

$ !${»}(*)) = {-iH(t) + z(t)L-J2g k L x (sm(n k t) 



fc>0 







+ caa(n h t)-J^+)}\ij>{ qk }(t)), 



dYr 



'9X+ 



(4.27) 



1. The Markov Limit 

Finally in this subsection we determine the Markov 
limit of this equation. Applying the results at the end of 
Sec. IV A, we get 

dt\A(t)) = (-»#(*) + Lz(t) - lL x £j \Mt)), (4.32) 

as O z (t,t) = L. This is in Stratonovich from. Wc 
transform this to the Ito form by using the method in 
Sec. MB 1. In this case the Ito correction is 



/./,-. 



(4.33) 



and the Ito SSE is 



d\Mt)) =dt[-^^+Lz(t)-^DL)\i, z (t)), (4.34) 
which is the general linear homodyne SSE [26, 33] as 



C. The Actual Stochastic Schrodinger Equation for 
the Quadrature Unraveling 



which is a linear equation solely in terms of the parame- 
ters {X+} and {Y~}. 

As in the coherent case, to make progress towards a 
genuine SSE we wish to replace the partial derivatives by 
a functional derivative with respect to the noise function. 
To do this we note that. 







by;. 



6 dzjs) 
5z(s) dX+ 

1 S dzjs) 
o Sz(s) dYr 



ds, 



ds . 



(4.28) 
(4.29) 



Thus we obtain 

d t \Mt)) = {-iH(t)+z(t)L-L x J p(t-s) 

X JzTs) ds }^ z(t)) > (4 ' 30) 

where f3(t — s) is the memory function for the noise. As 
in the coherent state case, this enforces an initial vacuum 
state for the bath. The final step to obtaining the linear 
non-Markovian SSE with real noise is to assume that the 
functional derivative can be replaced by an operator as 
in Eq. (3.23). With this Ansatz the linear SSE becomes 

d t \Mt)) = (-iH(t)+z(t)L-L x f p(t -s) 



xd z (t, s )ds)\Mt)) 



(4.31) 



As in the coherent case, to find an actual SSE (gener- 
ating states with the actual probability) we need to find 
random variables with the actual probabilities of mea- 
surement results {q k }- To sort these out we use the Gir- 
sanov transform (2.21) to find a first-order partial differ- 
ential equation (PDE) for the probability, from which the 
characteristic equation generates the transformed vari- 
ables 



d t P({x+,Yr},t) = {(4> {qk} (t)\dS {qk} (t))+c. 



xA({X+,Yr}). (4.35) 
Using equations (4.27) allows us to write 

d t P({X+,Yr},t) = -Y^g k ^+(cos(n k t)(i> {qk} (t)\ 



x£ K |V> te} (i))A({X+ y fe -}), 
~ zZ 9k rr^ ( sin(ft fe t) (ijj {qk } (t)\ 

k>0 ° Y k 

xL x \i> {qk} (t))A({X+,Yr}) 



(4.36) 



This can be simplified to 
dtP({X+,Y k },t) = 

fc>0 k 



(Q k t)(L x ) t P({X+,Y k -},t)) 
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T,9k^(^m)(L x ) t P({X+ : Y-},t) 



k>o 01 k 



(4.37) 



With the above random variable equations for Xt(t) 
and Y k ~ (t) we can write the noise function for the actual 
probability as 



where (L x ) t is defined by equation (3.32). 
The characteristic equations are 



— X+ = g k cos(tt k t)(L x 



dt k 



g k sm.{£l k t){L x ) t . 



(4.38) 
(4.39) 



Integrating these differential equation from time to t 
we get 



z(t) = z A {t)+ / (L x ) t (3(t-s)ds, (4.43) 
Jo 



where Z\ (t) is the random variable with statistics deter- 
mined by the A({X£ , Y k ~}) distribution. That is, the 
correlations of z\(i) are those of z(t) in Eq. (4.18). 



X+(t) = X+(0)+ / g k cos(n k s){L x ) s ds, (4.40) 
Jo 

**"(*) = Y k~(°)+ f 9kM^k-s)(L x ) s ds. (4.41) 
Jo 

The distribution for X k (0) and Y k ~ (0) is due to the quan- 
tum initial conditions. As before, the use of the func- 
tional derivative in Eq. (4.30) implies that the initial bath 
state is a vacuum state. Thus, the randomness in X k (0) 
and Yr(0) is that of the ostensible distribution: 



p«x+ y-},o) = A({x+ y-}) 



a Z—tk>0^ k ' k ' 



r«/2 



(4.42) 



Now we have the correct noise function we can cal- 
culate the actual SSE. As in the coherent case we need 
dt\^zit))i and for this case equation (2.29) will be 



d t \i> z {t)) = {-iH(t) + Lz{t)-{L x ~(L x ) t ) 

5 

] Jzis) 



o 



Following the same procedure as in the coherent case we 
obtain 



-iH{t) + {L-{L) t )z{t))\^ z {t)) 
1 



1 (L x -(L x ) t ) [ p( t -s) T ^ds\Mt)) 



5z(s) 



((L x - (L x ) t ) / p(t-s) J - r -rds\^ z (t))\Mt))- 
H{ qk }(t)\ Jo Sz{s) 



(4.45) 



Again this is not a SSE until we make the Ansatz defined 
in Eq. (3.23), which gives 



-iH(t) + (L-(L) t )z(t) 
[L x -{L x ) t ) f f3(t-s)6 z (t,s)ds 



+ ((L, - (L x ) t ) j p(t-s)O z (t,s)ds 



x )J\Mt))- 



(4.46) 



This is the actual SSE for real-valued noise. All of the 
comments regarding the interpretation of the correspond- 
ing complex-valued noise SSE (3.44) carry over to this 
case. 



1. The Markov Limit 



Taking the Markov limit of the actual SSE results in a 
noise function of the form 



z(t) = z A (t) + l(L x ) u 



where z^(t) = ^y£,(t). The actual SSE becomes 



(4.47) 



-zH(t) + (L-(L) t )(z(t) + j(L x ) t ) 



-±{L x L-{L x L) t ))W z {t)), 



(4.48) 



This is in Stratonovich form, to compare it to the equiv- 
alent homodync SSE we need to convert it to Ito form. 
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The Ito correction term for this equation is 

* y ( b ° bj + bt JL bj ) = *2 (££ - 2L(L) t 

-(L x L) t + (L x ) t (L) t + (L) t (L) t ^ z (t)). 

(4.49) 

Adding this to the Stratonovich SSE we get the following 
Ito SSE, 

d\i> z {t)) = dt{-iH(t) + (L-{L) t )z(t)-^dt(LiL 

-t(tf) t + L(L) t -{L) t {l) t )}\il>&)). 

(4.50) 

This is the same as the homodyne SSE presented in Ref. 
[25, 36] when we substitute in Eq. (4.47) for z(t). As in 
the coherent case there will be a difference between z(t) 
and the homodyne current, which from reference [25] is 

1(f) = \fl £,(t) + l{L x )t- This difference again comes 
down to the fact the in the quantum trajectory theory 
the measurement occurs a time dt later. 



V. A SIMPLE SYSTEM 

In this section we apply the above theory to a very 
simple non-Markovian system: a TLA coupled linearly 
and with the same strength to two single mode fields 
(labeled by k = ±1) that are detuned from u>o by ±A 
respectively. Without loss of generality, we can take the 
coupling strength g\ = g to be real. Then the memory 
function becomes 

a(t - s) = 2g 2 cos(A(t~ s)). (5.1) 

Note that this memory never decays, indicating that the 
dynamics of the atom is extremely non-Markovian. This 
is different from all cases considered by DSG, where the 
memory was taken to decay exponentially It is thus 
interesting to see how the formalism copes with this ex- 
treme case. At the same time, the simplicity of the bath 
(two modes) means that an exact numerical solution for 
Prcd(i) is relatively easy to find. This allows verification 
of the validity of the SSEs in reproducing p le df) by en- 
semble average, for both the linear and actual (nonlinear) 
cases. 

We would also like to see the different individual be- 
haviour of the trajectories corresponding to two different 
measurements (coherent state and quadrature measure- 
ments). This is readily apparent in this system for the 
initial condition |^(0)) = |e), where |e) and \b) are the 
excited and ground state of the TLA, respectively, so we 
choose this for all our simulations. 



A. Exact Solution 

To calculate the exact p ro d(t) we need to solve the 
Schrodingcr equation, which is displayed in Eq. (2.8). 
For this simple system we assume H = and 

V(t) =ge iA \a\a-a^)+ge- iAt (al 1 a-a 1 ^) (5.2) 

as S7i = A = — and g = g~\ = g\. Here the Lindblad 
operator L = a = \b)(e\. Since initially the field is in 
the vacuum state (|0i) (g> |0_i)) then the only non-zero 
complex amplitudes in are 

= ci(t)|600) + ca(t)|e00) + c 3 (*)|&01) + c 4 (i)|M0), 

(5.3) 

where |600) is short hand for \b) <g) |0i) ® |0_i) etc. Ap- 
plying the above Hamiltonian to this state we get the 
following four differential equations for the complex am- 
plitudes, 



ci(t) = 0, (5.4) 

c 2 (t) = - C3 (%e iAt - C4 (t) 5e - jAt , (5.5) 

cs(t) = c 2 {t)ge- lA \ (5.6) 

U(t) = c 2 {t)ge lA \ (5.7) 



which can be solved numerically For the initial state 
|e00), 02(0) = 1 and the rest are zero. Once we have the 
amplitudes for all time we know and by Eq. (1.1) 

we can then calculate p re d(t). For the TLA it is conve- 
nient to define the reduced state in terms of a pseudo-spin 
vector (x, y, z) by 

Pred{t) = -[/+ x{t)a x + y{t)a y + z{t)a z ], (5.8) 

where x(t), y(t) and z(t) are real parameters which equal 
the expected value of the corresponding spin matrix. 
These can be found from the above complex amplitudes 
by 

/ = | C iW| 2 + |c2(t)| 2 + | C3 W| 2 + |c4(t)| 2 , (5.9) 
x(t) = c 2 (t)c* 1 (t) + c* 2 (t)c 1 (t), (5.10) 
y(t) = -ic 2 (t)cUt)+ic* 2 (t) Cl (t), (5.11) 
z(t) = |c 2 (0| 2 -|c 1 (t)| 2 -|c 3 (0| 2 -|c4W| 2 -(5.12) 

To graphically illustrate the reduced state we numerically 
calculated the above real parameters for A = 2g. The 
results arc shown in Fig. 1 as a solid line. 



B. Coherent Unraveling 

For the simple system the memory function, Eq. (3.7) 
is given by Eq. (5.1), and and the noise operator for the 
coherent unraveling is, 

z(t) = 5 aie- lAt + gd^e lAt . (5.13) 
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FIG. 1: This figure depicts the reduced state calculated by 
three different methods; the exact solution (solid line), the 
ensemble average of 1000 SSEs for both the linear (dotted) 
and actual (dashed) SSE for the coherent unraveling. In this 
figure all calculations were done using a simple Euler method 
with a step size of dt — 0.0001, a detuning of A = 2g and 
initial system state of the form |^(0)) = |e). 



The linear SSE was obtained when we assumed an osten- 
sible probability A(ai,a_i) equal to the vacuum distri- 
bution 



A(ai,o_i) =7r - 2 e -l ai l 2 -l a - 1 l 2 



(5.14) 



With this probability distribution, we can write the noise 
function as a random variable equation of the form, 



z(t)=ga 1 er^ t +ga- 1 e^ t , 



(5.15) 



where a\ and a_i are complex GRVs of mean and vari- 
ance 1. 

Applying the simple systems dynamics to equation 
(3.22) we obtain 

dS z (t)) = (*•(*)* - & J «(* - s )j^) ds ) ^ z{t)) - 

(5.16) 

In Sec. Ill B we made the general Ansatz described by 
Eq. (3.23). For this simple system the specific Ansatz we 
will use is 

!&(*)> = /(*,*)*!&(*)>• (5-17) 



5z*{s) 



To work out the functions f(t, s) we use the following 
consistency condition [21], 

S !&(*)> = Jx7I^I&(*)>. ( 5 - 18 ) 



Sz*(s)dt^ zy " dtSz*{s) 
This gives 

d t f(t, s)&\^ z (t)) = /(*, s)F(t)&$ x {t)), (5.19) 



where 

F(t)= [ a(t-s)f(t,s)ds. (5.20) 
Jo 

This allows us to write the linear SSE for the coherent 
unraveling as 



d t \Mt)) = z'(t)*-**&F(t) \Mt)) 



(5.21) 



This is simple to to solve numerically provide we have a 
solution for F(t). 

The best way to calculate F(t) is to split it into to two 
terms, F(t) = Fi(t) + F_i(t) where, 



\g\ 2 e- iA ^f(t,s)ds = F*_ 1 (t). (5.22) 



Differentiating the above equations for Fi(t) and F_i(t) 
and using Eq. (5.19) and the fact that f(t,t) = 1 yields 

d t F 1 (t) = \g\ 2 -iAF 1 (t)+F 1 (t)F(t), (5.23) 
d t F^(t) = \g\ 2 +iAF_i(t) + F_i(t)F(t), (5.24) 

which can be solved numerically The initial conditions 
are F(0) = Fi(0) = F_i(0) = 0. Writing \i> z (t)) = 
C e (t)\e) + Cb(i)\b) gives us the following two differential 
equations, 



d t C e (t) = -C e (t)F(t), 
d t C b (t) = z*{t)C e {t). 



(5.25) 
(5.26) 



For an excited-state initial condition these equation can 
be solve numerically. Note that these solutions will not 
remain normalized, and the norm of most of them be- 
comes very small. This reflects the fact that a typical 
individual solution of this SSE does not correspond to 
a typical measurement result. Nevertheless, the ensem- 
ble average of the unnormalized states is p Te d{t)- To show 
this we simulated 1000 SSE for different z(t). The results 
of this simulation are shown in Fig. 1 as a dotted line, 
where the agreement with the exact solution is good. 

The actual SSE for coherent unraveling is found by 
applying the above results to Eq. (3.44). Doing this we 
obtain 

d t \i>S)) = {(v-(<?) t )z*(t)-{^ -(^) t )aF{t) 

+ ((^-(^) t )a)F(t)}\Mt))- (5-27) 
The noise, z*(t) in this equation is given by, 

z*{t) = z* A {t) + f a*(t - s)(^) s ds, (5.28) 
Jo 

where z\ (t) is the noise function used in the linear case. 
With this SSE the two differential equations for the com- 
plex amplitudes become 

d t c e (t) = -c!(t)c;(t) z *(t) + F(t)c e (t)(-i 

+ \C e (t)\ 2 -\C e (t)\ 2 \C b (t)\ 2 ), (5.29) 
d t C b (t) = C e (t)(l-\C b (t)\ 2 )z*(t)+F(t)C b (t)\C e (t)\ 2 
x(2~\C b (t)\ 2 ). (5.30) 
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FIG. 2: This figure shows a typical trajectory generated by 
the actual SSE for both the coherent (solid) and quadrature 
(dotted) unraveling. These were all done with the parameters 
defined in Fig. 1. 

The solution to these equations is an actual state, in 
the sense that it is normalized, and generated with the 
actual probabilities. Thus a typical trajectory docs give, 
at any time t, a typical state that corresponds to an ob- 
server measuring it at that time in the coherent basis. It 
is thus worth examining a typical trajectory, which we 
have plotted in Fig. 2 (the solid line). The normaliza- 
tion of the state is shown to remain equal to one, within 
the error introduced by the integration algorithm. To 
show that the ensemble average of these trajectories is 
the reduced state, an ensemble average of 1000 SSE was 
simulated and the results are depicted in Fig. 1 (dashed 
line). We see that the actual case is closer to the p re d(t) 
then the linear case. This is expected as in general the 
linear SSE converges slower than the actual SSE, as most 
of the states generated from the linear SSE have virtually 
no contribution to the mean. 



C. Quadrature Unraveling 

If we apply the theory for the quadrature unraveling to 
this simple system, the quadrature noise operator, equa- 
tion (4.4) becomes 

z(t) = 2g{X+ cos(At)) + Yf sin(At)}, (5.31) 

and the quadrature noise function is 

z{t) = 2g{X+ cos(At) + Ff sin(At)}, (5.32) 

which is real. If we choose the ostensible probability to 
equal the vacuum probability, then 

A(X+, Ff ) = 7r- 1 e- x * a - Y i~ 2 . (5.33) 



Thus for the linear case X^ and Y[~ are GRVs of mean 
zero and variance 1/2. 

For this simple system the quadrature linear SSE, 
Eq. (4.30), becomes 

d t \Mt)) = (z(t)a-a x J p( t -s)-^-rds)$ z (t)) 

(5.34) 

As for the coherent case we can make an Ansatz for the 
functional derivative. We again choose Eq. (5.17). This 
allows us to write the quadrature linear SSE as 

dt\$,(t)) = (z{t)& - a x aF{t)) |&(t)), (5.35) 

where F(t) is given by 

F(t) = f (3(t-s)f(t,s)ds, (5.36) 
Jo 

and (i(t - s) = 2\g\ 2 cos(A(t - s)). 

It turns out for this simple system F(t) is the same for 
both the coherent and quadrature unraveling, because 
a(t — s) = j3(t — s). Knowing F(t), we get the following 
two differential equations for the state: 

d t C e (t) = -C e (t)F{t), (5.37) 
d t C b (t) = z(t)C e (t). (5.38) 

These are the same as for the coherent case, except that 
z(t) is generated differently. To show that the ensemble 
average of the solution to the linear SSE for the quadra- 
ture unraveling converges to /3 re d(i), 1000 trajectories for 
different z(t) where simulated. The results of these sim- 
ulations are shown in Fig. 3 as a dotted line, where it 
is seen that the ensemble average of the linear SSE does 
reproduce the exact solution for p Te d(t) with little error. 

The actual SSE for quadrature unraveling is found by 
applying the above results to Eq. (4.46), 

dMM) = {(a - (a) t )z(t) - (a x - (a x ) t )aF(t) 

+ - (ff x )t)ff) F(t)}|V>«(t)>. (5-39) 
The noise, z(t) in this equation is given by, 

z(t) = z A (t)+ f /3(t-s){a x } s ds, (5.40) 
Jo 

where z\(t) is the noise function used in the linear case. 
With this SSE the two differential equations for the com- 
plex amplitudes become, 

d t C e (t) - F(t)C e (t)(-l + \C e (t)f - \C e (t)\ 2 \C b {t)\ 2 ) 
-F{t)Cl{t)Cf{t) - Cl{t)Ct{t)z{t), (5.41) 

d t C b (t) = F(t)C b (t)\C e (t)\ 2 (2-\C b (t)\ 2 ) 

+F(t)C* b (t)C 2 e (t)(l - \C b {t)\ 2 ) + C e (t) 
x(l-\C b (t)\ 2 )z(t). (5.42) 
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FIG. 3: This figure depicts the reduced state calculated by 
three different methods; the exact solution (solid line), the 
ensemble average of 1000 SSEs for both the linear (dotted) 
and actual (dashed) SSE for the quadrature unraveling. These 
where all done with the parameters defined in Fig. 1. 



A typical trajectory from the quadrature SSE is illus- 
trated in Fig. 2 (the dotted line). Note the feature that 
clearly distinguishes it from the coherent trajectory: y is 
always zero. To show that the solution of the actual SSE 
reproduces the reduced state on average, an ensemble of 
1000 actual SSEs was simulated and the results are de- 
picted in Fig. 3 (dashed line). We see that it reproduces 
the exact solution, again with less error than that from 
the linear SSE. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we have explored non-Markovian stochas- 
tic Schrodinger equations by furthering the work of Diosi, 
Strunz and Gisin [16, 19, 20, 21, 22]. Specifically, we have 
interpreted their results in the framework of quantum 
measurement theory. Their SSEs arise as a special case 
when the measurement basis of the bath is the coherent 
states, so we label it the coherent unraveling. The benefit 
of using the measurement interpretation is two-fold. 

First, it allows us a better understanding of the inter- 
pretation of non-Markovian SSEs. The state at any time 
t generated by the SSE can be interpreted as a condi- 
tioned system state, given a particular result from a par- 
ticular measurement on the bath. However, the measure- 
ments at different times are incompatible, so the linking 
together of different states over time is, we have argued, 
a convenient fiction. Thus the trajectory generated by 
a non-Markovian SSE does not have the same physical 
status as that generated by a Markovian SSE, where the 
measurements at different times arc compatible and the 
states at different times can represent a single evolving 
system. 



Second, it allows us to generate other sorts of SSEs 
corresponding to different sorts of measurements on the 
bath (unravelings). In this paper we presented a second 
unraveling, based on measuring certain quadrature oper- 
ators on the bath. This gives rise to an SSE only under 
certain assumptions to do with the bath frequencies and 
couplings. The resultant SSE contains real-valued noise, 
as opposed to the complex noise in the SSE of DSG. 
The ability to construct a non-Markovian SSE with real- 
valued noise is contrary to the expectation expressed by 
DSG in [21]. 

We have also shown in this paper that the Markov 
limit of the quadrature and coherent unravelings are ho- 
modyne and heterodyne detection respectively. As noted 
above, in this Markov limit the SSE generates a true 
quantum trajectory for a conditioned system state over 
time. It is interesting that this arises smoothly as the 
limit of a non-Markovian SSE that does not have this 
interpretation. However, as we have shown, one has to 
be very careful with the definition of the time of mea- 
surement in order to reconcile this limit with the usual 
quantum trajectory theory. 

To illustrate our general theory we have applied it to a 
simple system: a TLA coupled linearly to just two single- 
mode fields detuned from the atom by ±A. This is an 
extremely non-Markovian problem with no finite memory 
time, unlike the previous examples considered by DSG. 
Nevertheless the theory is able to describe the evolution 
of the atom by an SSE. In Fig. 2 we displayed typical 
non-Markovian SSEs for both the quadrature and coher- 
ent unraveling, and in Figs. 1 and 3 we showed that on 
average both SSEs do generate the exact reduced state. 

In conclusion this paper has presented a significant 
generalization of the DSG approach to non-Markovian 
SSE. However, there is still a lot of questions to be an- 
swered. 

First, is it possible within this framework to derive 
other classes of non-Markovian SSEs? In particular, is 
it possible to describe an unraveling based on discrete 
measurement on the bath, say the in number-state basis? 

Second, is there a physical system where our theory 
could be naturally applied? That is, is there a physical 
system where the bath could be measured in a suitable 
basis at an arbitrary time so as to produce a pure condi- 
tioned system state? 

Third, what conditions are necessary for one to be 
able to find a suitable Ansatz for replacing the functional 
derivative with an operator? As we have argued, this is 
necessary to create a genuine SSE. Yu, Diosi, Gisin and 
Strunz have given a general procedure for finding this op- 
erator, but only when the system dynamics are weakly 
non-Markovian (the so-called 'post-Markovian' approx- 
imation) [39, 40]. We suspect that the conditions for 
finding an exact Ansatz depend upon both the nature of 
the system and its coupling to the bath. 

Fourth, can the techniques of non-Markovian SSEs be 
applied as a numerical tool for studying real systems? 
We have in mind potentially strongly non-Markovian sys- 
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terns such as an atom laser [41] or photon emission in a 
photonic band-gap material [42, 43]. 

Fifth, and last, is there an alternative framework to 
standard quantum measurement theory in which there 
is a physical interpretation for a trajectory generated by 



a non-Markovian SSE? That is, can the states at differ- 
ent times in a single trajectory generated by the SSE be 
interpreted as pertaining to a single system in some non- 
standard approach to quantum measurements? This is a 
very open question. 
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